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ABSTRACT 

SPACE  AND  ENERGY  DEPENDENT  NEUTRON 
SLOWING  DOWN  TIME  IN  GRAPHITE 


Submitted  to  the  Department  of  Nuclear  Engineering  on 
September  26,  1967  in  partial  fulfillment  of  the  requirements 
for  the  degree  of  Master  of  Science. 

A  4  foot  cube,  graphit^  moderating  assembly  was  con- 
structed at  the  M.I.T.  Rockefeller  Van  de  Graaff  accelerator. 
Using  the  pulsed  neutron  experimental  installation,  a  neu-   . 
tron  burst  of  200  nsec  width  was  injected  into  the  center 
of  the  assembly.   Prompt  gamma  radiation  from  a  neutron 
resonance  capture  by  an  indicator  foil  in  the  moderator 
was  measured,  with  a  3  inch  x  3  inch  sodium  iodide  scin- 
tillation detector,  as  a  function  of  time  after  the  neu-     ; 
tron  pulse.   The  energy  dependence  of  slowing  down  time 
was  experimentally  measured  by  using  resonance  absorber     ; 
foils  of  cadmium,  indium,  and  gold  with  isolated  reson- 
ances of  0.2,  1.46  and  5*0  eV  respectively.   By  vary- 
ing the  distance  of  the  foil  from  the  source,  the  spatial 
dependence  of  slowing  down  time  was  also  investigated.' 
Comparisons  are  made  with  slowing  down  theory.   The  experi-  : 
mental  results,  together  with  their  estimated  uncertainties,' 
are  tabulated  below: 


tm  (usee) 
r(cm)  0.2  1.46  3.0 

12.5      47s.  50  +  3.00     19.75  +  .50  11.25  +  .25 

18.0      47.25  +  2.25     19.37  +  .25  10.10  +  .25 

28.0      47.25  +  2.25     19.00  +  .50  10.87  +  .25 

37.5  -  20.00  +  .50  11.13  +  .37 

44.0  -  -  12.50  +  .25 

44.5      49.75  +  2.25     21.62  +  .50 

Thesis  Supervisor:   Franklyn  M.  Clikeman 

Title:   Assistant  Professor  of  Nuclear ' Engineering 
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CHAPTER  1 
Introduction 

In  his.M.I.T.  Ph.D.  thesis,  Adamant iades  (l)  demon- 
strated  that  the  pulsed  neutron  technique  using  a  resonance 
absorber  as  a  point  indicator  could  be  used  to  measure  sys- 
tem transient  effects  such'  as  the  slowing  down  time.   He 
designed- and  installed  an  experimental  installation  for 
pulsed  .neutron  experiments,  and  made  some  measurements  in 
water.   Vaughan  (2)  refined  Adamant aides'  system,  and  measured 
the  spatial  dependence  of  slowing  down  time  in  water  for 
three  different  resonance  absorbers.   Moller  has  carried 
out  extensive  experiments  in  the  field  of  neutron  slowing 
down  for  both  light  (5,4)  and  heavy  (5)  water.   However, 
there  does  not  appear  to  be  any  published  measurements  of 
the  slowing  down  time  for  heavier  moderators.   Therefore, 
in  this  present  experiment,  the  pulsed  technique  has  been 
expanded  to  investigate  slowing  down  time  in  graphite. 

In  the  pulsed  neutron  method,  a  short  monoenergetic 
burst  of  neutrons  is  injected  into  the  moderating  medium 
at  a  point.   The  neutrons,  originally  high  energy,  gradually 
slow  down  by  elastic  collisions  with  the  moderating  medium. 
If  a  neutron  indicator  is  incorporated  into  the  moderating 
medium,  it  is  possible  to  measure  the  time  difference  between 
pulse  injection  and  neutron  arrival  at  the  indicator.   A 


aac 


i 
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time  distribution  is  formed  from  many  time  differences, 

measured  under  the  exact  same  conditions. 
/ 

Since  most  materials  emit  prompt  gamma  rays  upon  the 
capture  of  a  neutron,  detection  of  these  gamma  rays  pro- 
vides one  method  by  v/hich  the  time  behavior  of  neutron 
captures  may  be  followed  after  the  pulse.   If  an  absorber 
with  a  strong,  isolated  neutron  capture  resonance  is  incor- 
porated into  the  moderating  system  as  a  point  indicator, 
then  the-  reaction  rate  of  the  neutron  flux  with  the  absorp- 
tion resonance  is  proportional  to  the  time  dependent  flux 
at  that  point,  for  neutrons  with  energies  equal  to  the 
resonance  energy.   That  is,  the  neutron  density  (or  flux) 
N(r,E   t) 

N  -  neutron  density 

r  -  distance  from  source 

E  -  energy  of  absorber  resonance 

t  -  time 
can  be  measured.   Thus,  variation  of  the  position  of  the 
absorber  and  the  use  of  different  resonance  absorbers  should 
yield  the  time,  space,  and  energy  variation  of  the  neutron 
density  following  the  pulsing  of  the  assembly. 

By  using  the  M.I.T.  Rockefeller  Van  de  Graaff  accelera- 
tor as  a  source  of  high  energy  protons,  high  energy  neutrons 

7        7 
were  obtained  by  the  Li' (p,n)Be^  nuclear  reaction.   The  lithium 

target  was  located  at  the  center  of  a  large  assembly  of 
graphite,  thus  approximating  a  point  source  in  an  infi- 
nite medium.   The  accelerator  beam  was  pulsed  by  an  external.  _'. 
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pulsing  system  after  the  protons  had  been  accelerated.   The 
pulse  width  used  in  these  experiments  v/as  200  nsec.   Enough 
time  was . allowed:  to  elapse  between  pulses  so  that  all  tran- 
sient neutron  effects  from  one  pulse  would  die  out  before 
the  next  pulse  was  injected.   Absorber  foils  of  cadmium, 
indium,  and  gold,  with  isolated  resonances  located  at  0.2  eV, 
1.46  eV,  and  5«0  eV  respectively,  were  used  to  sample  the 
neutron  flux  at  these  energies.   A  detecting  and  analyzing 
system  was  assembled  that  was  capable  of  indicating  the 
time  difference  between  pulse  injection  and  the  resonance 
capture  of  a  neutron.   The  most  probable  time  difference 
from  injection  to  resonance  capture  was  assumed  to  be  the 
neutron  slowing  down  time  for  that  energy  at  that  position. 

The  material  in  this  thesis  is  organized  as;  follows: 
Chapter  2  is  concerned  with  slowing  dov/n  theory}  the 
experimental  facilities  are  described  in  Chapter  3|  the 
performance  of  the  experiments  is  described  in  Chapter' 4-* 
and  the  results  are  presented  and  discussed. 
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.  CHAPTER  2  i 

Theory 

|  '         I 

■2.1   General  ^ 

The  analytical  treatment  of  pulsed  neutron  assem- 
blies requires  a  mathematical  solution  of  the  general 
transport  equation.   If  an' expression  representing  the  neu- 
tron density  as  a  function  of  time,  energy,  and  space,  and 
satisfying  the  Boltzmann  transport  equation  with  its  source 
and  "boundary  conditions  can  be  found,  the  behavior  of  neu- 
trons in  the  "system  can  be  completely  described. 

Unfortunately,  no  exact  solution  of  the  Boltzmann 
equation,  in  all  its  generality,  has  ever  been  found^ 
but  approximate  solutions  are  available.   Of  these  solutions, 
',  only-. a  few  may  be  verified  because  not  all  quantities  can 
be  measured  experimentally. 

The  quantity  most  readily  measured  experimentally 

in  the  slowing  down  region  is  the  neutron  slowing  down 

time  (t  ).   This  is  the  time  at  v/hich  the  neutron  density 
m 

at  some  energy  (E)  is  maximum  and  can  be  found  by  setting 
the  derivative  of  N(r,E,t) 

N  -  neutron  density 

r-  distance  from  source 

E  -  neutron  energy 

t  -  time  •  ■  • 


. 
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with  respect  to  time  equal  to  zero. 

It  is  also  possible  to  measure  the  mean  slowing  down 
time  (<t>)  which  is  also  defined  at  the  first  time  moment. 
The  time  moments  are  integral  quantities  defined  as: 

^(r.v^-J^t11  N(v,t)  dt/f^N(v,t)  dt     (2.1.1) 

lT(v,t)  -neutron  density  as  a  function  of  space 
and  energy  (velocity) 

t  -  time 

Both  the  time  to  maximum  or  the  slov/ing  down  time, 
and  mean  slov/ing  down  time  may  "be  measured  experimentally  as 
a  function  of  space  and  neutron  energy. 

In  this  experiment,  the  source  may  be  considered 
isotropic  and  a  delta  function  in  time,  energy,  and  space. 
Although  the  neutron  source  energy  does  vary  with  the  angle 
in  the  laboratory  system,  this  source  term  can  be  neglected 
because  it  appears  in  the  denominator  of  the  expressions 
for  the  time  distributions,  and  it  is  much  greater  than  the 
energy  at  which  the  time  distribution  is  measured.   For 
example,  one  expression  for  the  neutron  slowing  down  flux 
in  hydrogen  is: 

cp(v,t)  =  S.(vt/ys)2  • exp(-vt/Ys)•(l-v/vo+2Ys/vot)+ 
2S.vo/v.5(v-vo).exp(-vot/Ys)      (2.1.2) 

cp  -  neutron  flux 

v  -  neutron  velocity 

S  -  source  strength 
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Y   -  scattering  mean  free  path 

v  -  neutron  source  velocity 
o  J 

t  -  time 

However,  because  v  »v,  equation  (2.1.2)  simplifies  to: 
cp(v,t)  =  S.(vt/YJ2.  exp(-vt/yj      (2.1.3) 
Except  for  positions  closer  than  three  mean  free 

paths  away  from  the  physical  boundary  of  the  assembly,  the 

spatial  boundary  conditions  are  those  of  an  infinite  medium. 

For  positions  closer  than  three  mean  paths,  leakage  must 

be  taken  into  account. 

In  the  energy  region  lying  above  a  few  tenths  of  an 

eV,  the  nuclei  of  the  moderating  assembly  can  be  considered 

•  ■"!  i 

free  and  at  rest.   For  lower  energies,  however,  the  molecular 

I; 

structure  of  the  moderator  becomes  important. 

Any  absorption  which  occurs  in  the  moderating  assembly 

will  effect  the  time  distribution.   Those  neutrons  which 

take  longer  to  slow  down,  i.e.  have  more  collisions,  to 

energy  (v)  will  have  a  greater  probability  of  being  absorbed. 

Thus,  absorption  contributes  to  a  decrease  in  t  .   In 
'  .        ■  •  m 

graphite,  the  absorption  cross  section  is  very  small^  and 
absorption  effects  may  be  neglected. 

2.2   Space  Independent  Slowing  Down  Theory 

If,  in  addition  to  an  infinite  medium,  we  assume  a 
uniformly  distributed  source  instead  of  a  spatial  delta 
function  source,  the  spatial  dependence  is  removed  from  the 
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Boltzmann  Equation.   The  equation  reduces  to: 
.  l/v-.6(p(v,t)/6t  =  -[s  (v)  +  s(v)]cp(v,t)  + 

is         ci 

•  Jcp(v',t).e(v' — >v)  dv'+S-6(v-v  )-6(t)      (2.2.1) 

cp(v,t)  -  neutron  flux  as  a  function  of  space  and 

energy  (velocity) 

• ,'   e    macroscopic  scattering  cross  section 

e   -  macroscooic  absorption  cross  section 
a  " 

.  e   -  removal  cross  section 

s  -. source  strength  (n/cm  -sec) 

v  -  source  energy 

6(  )  -  Dirac  delta  function 

Unfortunately,  for  a  heavier  than  water  moderator, 
it  is  very  difficult  to  obtain  an  expression  for  the  flux, 
and  hence,  for  the  slowing  down  time.   Using  a  method 
introduced  by  Marshak  (1)  the  time  moments  can  be  obtained 
exactly.   A  few  of  the  moments,  in  terms  of  the  moderator 
mass  number,  obtained  by  Adamantiades  (2)  are  given  in 
Table  I,  of  his  thesis.   Once  the  time  moments  have  been 
determined,  trial  functions  can  be  found  to  fit  these 
moments. 

Marshak  and  vonDardel  (3)  suggest  using  the  function: 

<p(x)  =  A.x"1+2/(l~p2).exp(-x-b/x)    •  (2.2.2) 

A  -  normalization  constant 

x  =  vt/y  where  y   -  scattering  mean  free  path 

r  =  (M-m)/(M+m) 

M  -  mass  of  moderator 
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m  -  mass  of  neutron 

"b  -  a  parameter 
which  gives  the  right  asymptotic  dependence  for  the  higher 
positive  moments.   vonDardel  varies  the  parameter  (b)  with 
(r)  wh'ile  Marshak  treats  it  as  a  constant. 

For  heavier  elements,  the  time  dependent  neutron  den- 
sity distribution  approaches  a  narrow  Gaussian  in  shape. 
Thus,  for  heavier  elements,  t   should  approximately  equal 
<t>. 

Adamantiades  gives  a  complete  discussion  of  the  time 
moment  method  in  Chapter  II  of  his  thesis.   He  also  men- 
tions other  methods,  and  gives  a  number  of  references  deal- 
ing with  the  problem  of  a  heavier  than  water  moderator. 

t: 

None  of  these  references,  however,  specifically  give  results 
for  graphite  assemblies. 

2.3   Space  Dependent  Slowing  Down  Theory 

When  the  neutron  source  is  localized,  as  in  this 
experiment,  the  effects  of  diffusion  during  moderation  must 
be  taken  into  account.   The  problem  becomes  much  more  com- 
plicated since  the  transport  equation  is  now  a  function  of 
three,  rather  than  two  variables. 

As  in  the  space  independent  case,  Adamantiades  in  Chap- 
ter II  of  his  thesis  discusses  the  theoretical  work  which 
has  been  done  in  this  area.   This  concludes  a  solution  by 
Dyad'. kin  and  Batalina  (4)  ,  for  the  space-energy-time  neu- 
tron density  distribution  assuming  no  absorption,  constant 
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scattering  cross  section,  and  isotropic  scattering.   The 

expression  involves  three  factors  which  are  complicated^ 

and  require  coefficients  which  are  dependent  on  the  medium 

being  examined.   The  authors  also  give  a  simplified  expres- 

i   sicn  for  t  ,  which  is: 
■  m* 

tmax  =  2Ys/vn.(2-3u/YFT^r).    .  ' 

[l-2a/3.(l-puArrrTpl7")']"1     (2.3.1) 
i  >  ■ 

Y   -  scattering  mean  free  path 

v  -  neutron  velocity 

7i  -  average  logarithmic  energy  decrement 

u  -  neutron  lethargy 

a  -  parameter  which  depends  on  the  moderating  medium 

p  -  parameter  which  depends  on  the  moderating  medium 

r  -  distance  from  source 


Results  are  given  for  sand  and  water  mixtures.  / 

Claesson  (5)  has  also  determined  the  space  dependent 
neutron  density  distribution  for  both  a  point  and  plane 
source  assuming  zero  absorption,  isotropic  scattering, 
and  constant  mean  free  path.   Although  his  method  can  be 
extended  to  moderators  of  arbitrary  mass,  curves  are  only 
given  for  the  case  of  water  at  the  indium  resonance.   These 
curves  show  however,  that  past  10  cm  the  slowing  down  time 
is  strongly  dependent  on  neutron  energy  at  birth. 

Adamantiades  also  develops  a  method  of  calculating, 
time  moments  using  an  iterative  solution  of  the  stationary 
transport  equation.   The  nth  time  mpment  is  expressed 
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as : 

.  <tn>  =  cp(n)/cp(0)      .(2.5.2) 

HcpCo)  =  -S(r,v,U)      (2.5.3) 

H(p(n)+n/v.cp(n-1}  =  0     (2.3.4) 

H  =  [-u^-Et+/dv'/dvL'  es(v'— *v,v.'— »u)]      (2. $.5) 

r,v,ji  -  vectors  for  position,  velocity  and  solid 
angle,  respectively 

e.  -  total  macroscopic  absorption  cross  section. 
Any  calculational  method  suitable  to  solve  the  steady  state 
can  be  used  to  solve  for  the  flux  moments. 

David  Diamond,  a  graduate  student  at  M.I.T.,  has 
kindly  supplied  -theoretical  time  distributions  at  the  gold 
and  indium  resonances  for  a  number  of  spatial  points, 
using  the  diffusion  approximation.   A  distribution  for  each 
foil  at  a  distance  from  the  source  of  18  cm  is  shown  in 
Figure  2.1  (a)  and  (b) . 
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CHAPTER  3 

Experimental  Equipment 

3»1   General  Description 

The  equipment  necessary  for  measuring  neutron  transient 
effects  by  the  pulsed  source  method  can  be  grouped  into  six 
categories: 

1)  The  neutron  source. 

2)  A  method  for  pulsing  the  source. 

3)  The  moderating  assembly. 

4)  A  system  for  -detecting,  analyzing  and  storing  data. 

5)  A  system  to  monitor  the  experimental  equipment. 

6)  A  normalizing  system  to  insure  consistent  data. 

The  general  layout  of  the  experimental  system  showing  I 
the  major  components  is  given  in  Figure  3.1. 

The  neutron  source  consisted  of  a  target  of  natural 
lithium,  located  approximately  in  the  center  of  the  modera- 
ting assembly,  bombarded  by  2.4-  MeV  protons  to  produce  the 
Li'(p,n)Be'  reaction.   The  high  energy  protons  were  produced 
by  a  Van  de  Graaff  electrostatic  accelerator.   The  pulsing 
system  electrostatically  deflected  the  proton  beam  off  the 
lithium  target,  allowing  it  to  strike  the  target  for  a  short 
time. 

Thus,  the  pulse  of  fast  neutrons  produced  at  the  target 
was  ejected  into  the  center  of  the  moderating  assembly. 
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These  neutrons  were  then  slowed  down  by  elastic  collisions 
in  the  medium,  and  their  time  distribution  was  measured 
by  the  detecting  and  analyzing  system.   This  system  con- 
sisted of  the  conventional  detector-preamplifier-amplifier 
array,  plus  electronic  equipment  for  pulse  height  discri- 
mination and  tine  analysis  combined  with  a  multi-channel 
analyzer  to  store  the  data.   There  was  a  common  link  betv/een 
the  pulsing  and  the  analyzing  system  which  provided  a  zero 
time  reference. 

The  monitoring  system  was  used  to  control  the  perfor- 
mance of  the  pulsing  system,  and  thus  to  optimize  the  quality 
of  the  neutron  pulse.   The  normalizing  system  yielded  the 
time  integrated  neutron  flux  for  each  experiment  performed. 
This  information  was  necessary  to  properly  carry  out  back- 
ground subtraction. 

A  more  detailed  description  of  the  experimental  equip- 
ment will  be  given  in  the  sections  to  follow. 


3.2  The  Mechanical  Layout 

A  plan  view  of  the  mechanical  layout  of  the  system  is 
given  in  Figure  3*2.   Only  a  brief  description  of  the  sys- 
tem will  be  given  here  since  a  thorough,  detailed  descrip- 
tion of  the  system  will  be  given  by  Adamantiades  (1),  who 
actually  constructed  the  system. 

Since  the  experimental  apparatus,  was  positioned  on 
the  floor  below  the  accelerator,  the  protons  were  accelerated 
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downward  and  bent  into  the  horizontal  bean  tube  by  a  bend- 
ing and  'analyzing  magnet.   The  beam  tube  was  kept  under  a 
high  vacuum  by  the  diffusion  pump  of  the  accelerator,  and 
by  a  small  Vac  Ion  pump  located  close  to  the  target. 

The  deflecting  plates  delivered  a  high  voltage  pulse 
which  deflected  the  beam.  The  water-cooled  chopping  slits 
were  mounted  after  the  deflection  plates  for  the  purpose 
of  intercepting  the  beam  when  it  was  swept  off  the  target. 
The  distance  between  the  slits,  the  shape  of  the  high  vol- 
tage pulse,  and  the  deflecting  plates  configuration  deter- 
mined the  shape  of  the  pulse  reaching  the  target. 

A  remotely  controlled  pneumatic  shutter  valve  was 
interposed  between  the  chopping  slits  and  the  target. 
The  beam  could  thus  be  cut  off  from  the  target,  enabling 

work  to  be  done  around  the  assembly.   The  target  tube  itself 

i 

was  constructed  of  aluminum  and  was  23  inches  long,  3/8 

i 

inch  O.D.,  with  a  1/16  inch  wall  thickness.   The  target 

i 

assembly  was  located  at  the  end  of  the  target  tube. 

It  was  decided  to  cool  the  target  in  order  to  protect 

it  from  deterioration.   This  was  done  by  sealing  the  end  cap, 

i 
which  presses  the  target  disc  directly  against  the  O-ring 

vacuum  seal,  with  a  plate  into  which  short  lengths  of  tubing 

had  been  brazed.   (See  Figure  3. 3)   V/ater' forced  through 

the  cap  removed  the  heat  from  the  target  disc. 

?  j 

3.3  The  Neutron  Source 

- 

As  in  previous  experimental  measurements  of  neutron 
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slowing  down  time  (1,2),  the  Li'(p,n)Be'  reaction  was  used 
as  the  neutron  source  mainly  on  the  basis  of  its  high 
yield  and  relatively  good  homogeneity.   For  a  given. angle 
from  the  target  and  for  laboratory  proton  energy  up  to  2.34- 

MeV,  the  reaction  is  monoenergetic.   The  threshold  of  the 

7        7^f 

Li'(p,n)Be'  reaction  is  shown  in  Figure  4.6  of  Adamantiades' 

thesis.   Because  of  neutron  intensity  considerations,'  a 
proton  energy  of  .2.4  MeV  was  used  in  this  experiment. 

The  target  used  in  these  experiments  consisted  of  a 
layer  of  natural  lithium  0.8/cui2  thick  deposited  through 
evaporation  on  a  circular  tantalum  disc  0.055  inches  thick 
and  13/16  inch  in  diameter.   Although  the  thickness  of  the 
target  is  less  than  1/10  of  the  range  of  2.5  MeV  protons 
in  lithium,  some  protons  will  undergo  slowing  down  colli- 
sions before  undergoing  nuclear  reactions.   The  energy 
spread  of  the  actual  beam  of  protons  incident  on  the  tar- 
get, however,  is  very  small  in  the  Rockefeller  Van  de  Graaff 
accelerator. 

To  summarize,  there  were  three  major  sources  of  neutron 
energy  inhomogeneity  in  this  experiment: 

1)  Angular  dependence  of  neutron  energy. 

2)  Excited  product  nucleus  reaction. 

3)  Energy  spread  of  the  reacting  protons. 

As  mentioned  in  section  2.1,  this  energy  inhomogeneity 
could  be  tolerated  since  the  slowing  down  time  is  relatively 
insensitive  to  the  source  energy. 


c 
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3.4  The  Moderating;  Assembly 

The  moderating  assembly  (see  Figure  3«4)  consisted  of 
reactor-grade  graphite  stringers,  approximately  4.0  inches 
by  4.25  inches  and  4.0  feet  long,  stacked  so  that  alternate 
JLayers  were  at  right  angles  to  each  other.   The  overall 
dimensions  of  the  assembly  was  3.6  feet  by  4.0  feet  by 
4.0  feet.   The  entire  assembly  was  covered  with  aluminum 
foil  to  shield  the  monitoring  and  detection,  systems  from 
R.F.  pickup  coming  from  the  pulsing  system,  although  this 
pickup  was  never  entirely  eliminated. 

In  a  stringer  located  at  the  center  of  one  face  of 
the  pile,  a  hole  was  bored  so  that  the  target,  along  with  the 
cooling  water  tubes  and  the  cables  of  the  monitoring  sys- 
tem, could  be  inserted  into  the  center  of  the  assembly. 

"At  a  right  angle  to,  and  on  the  same  horizontal  plane 
as  the  beam  tube  was  an  open  channel,  4.0  inches  wide, 
4.25  inches  high,  and  19.3  inches  long,  which  led  from  the 
target  to  one  face  of  the  assembly.   It  was  through  this 
channel  that  the  detector,  foil,  and  graphite  spacers  were 
inserted  or  removed  (see  Section  4.2). 

3.5  Electronic  Systems 
3.3«1   Beam  Pulsing  System 

Although  the  Rockefeller  accelerator  has  a  terminal 
pulsing  capability,  its  single  sweeping  frequency  was  much 
too  high  for  this  experiment.   Therefore,  a  post-accelera- 
tion system  was  used  which  had  the  advantage  of  a  variable 
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pulse  width  and  repetition,  rate  (see  Figure  3.5).   A 
detailed  discussion  of  the  system  is  given  in  Reference  (1). 
A  repetition  rate  of  2kc/sec  was  used,  allowing  500  usee 
between  pulses,  which  was  enough  time  to  allow  all  tran- 
sient neutron  levels  to  die  out.   A  pulse  width  of  200  nsec 
was  used  in  this  experiment. 

The  pulsing  system  was  triggered  by  the  general  Radio 
wave  generator  (type  no'.  1217-A)  ,  which  was  also  used  to 
start  the  timing  circuits  of  the  analyzing  system  (see 
Section. 3. 5«3) •   The  output  pulse  of  the  G.R.  pulser  ' 
triggered  an  E-H  Research  Laboratories  Model  No.  130  pulser. 

The  variable  width  output  pulse  of  the  E-H  pulser  was 
the  input  for  the  driver  circuit  which  is  essentially  a 
pulse  amplifier.   An  extensive  discussion  of  the  driver 
circuit  is  given  by  A.E.  Waltar'  (3)»  who  designed  it. 
Vaughan  modified  the  circuit  somewhat  to  give  a  narrower 
pulse.   In  this  experiment  the  driver  circuit  delivered  a 
2.9  KeV  pulse  to  the  deflection  plates. 

3>5»2   The  Detection  System 

The  object  of  these  experiments  was  to  detect  neutron 
capture  gamma  rays  as  a  function  of  time,  so  an  efficient 
gamma-ray  detector  was  required.   A  Harshaw  Integral  Line 
Scintillator  Assembly,  consisting  of  a  3  inch  x  3  inch 
Nal(Tl)  crystal  attached  to  a  6363  DuMont  photomutiplier 
tube,  was  used.   The  crystal  and  tube  are  hermetrically 
sealed  within  an  aluminum  housing,  making  the  detector 
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a  sin:    plug- in  unit.   A  thermal  neutron  shield,  consisting 
of  3  mm.  thick  boron  carbide  enclosed  "by  aluminum,  was  pro- 
vided to'  keep  activation  of  the  sodium  at  a  minimum. 

The  solid  state  preamplifier,  v/hich  plugged  into  the 

-end  of  the  detector  assembly,  was  designed  and  built  in  a 
cylindrical  can  having  the  same  diameter  as  the  detector. 
Vaughan  gives  a  schematic  diagram  of  the  preamplifier  in 
Figure  4.6  of  his  thesis. 

The  high  voltage  for  the  phototube  was  provided  by 
a  2  KV  power  supply,  while  the  preamplifier  was  powered  by 
a  small  solid  state  power  supply.   The  entire  detector 

'was  wrapped  in -aluminum  foil,  and  all  leads  to  and  from 
the  detector  were  enclosed  within  copper  braid  to  eliminate 
R.F.  pickup.  . 

$.5.3  The  Analyzing;  System 

The  main  features  of  the  time  analyzing  system  were: 

1)  A  start  pulse  providing  a  time  reference,  derived 
from  the  beam  pulsing,  network. 

2)  A  stop  pulse  indicating  the  time  at  which  an  event 
of  interest  took  place,  derived  from  the  detector. 

3)  A  timing  unit,  used  in  conjunction  with  a  multi- 
channel analyzer,  which  produced  a  count  in  the  channel 
corresponding  to  the  time  interval  between  the  arrival  of 
the  start  and  stop  pulses.   Conventional,  commercially 
available  equipment  was  used  in  the  entire  analyzing  system 
(see  Figure  3.6) . 
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The  start  pulse  was  obtained  from  the  General  Radio 
trigger  (see  Section  3. 5.1).   This  pulse  was  fed  through 
a  General  Radio  1392-A  Time  Delay  Generator  which  provided 
a  stable,  variable  tine  delay.   Unfortunately,  the  .output 
pulse  of  the  G.R.  delay  was  too  narrow  to  trigger  the  tim- 
ing unit',  and  an  integrating  circuit  was  employed  to 
stretch  the  pulse. 

The  detector  output  signal  was  fed  to  two  Ortec  4-10 
Strobbed  Discriminator,  with  the  strobe  input  coming  from 
an  Ortec  40?  Zero  Cross-over  Pickoff  Unit.   The  cross- 
over pickoff  v/as  fed  by  the  other  Ortec  4-10,  whose  time  con- 
stants had  been  optimized  for  a  stable  cross-over  point. 

i 

A  positive  logic  signal  was  provided  by  the  discriminator    ! 

i 

when  the  incoming  detector  signal  met  minimum  energy  require- 

i 
ments.   The  strobe  assured  a  constant  time  relation  between 

the  detector  input  and  discriminator  output.   The  positive    j 

logic  pulse  was  attenuated  by  a  resistor  network  before  it 

i 

was  fed  into  the  timing  unit. 

The  timing  unit  consisted  of  a  TMC  Model  211  Time  of 
Flight  Logic  Unit  which  plugged  into  the  TMC  Model  CN-110 
256  Channel  Pulse  Analyzer.   The  211  unit  converted  the  time 
difference  between  the  trigger  input  and  a  detector  signal 
to  an  appropriate  channel  number,  depending  on  the  scale 
factor  which  could  be  varied  from  0.25  to  64.  p.sec/channel. 
This  information  was  then  stored  in  the  multi-channel  analy- 
zer. 

The  TMC  211  Unit  could  handle  more  than  one  stop  pulse 
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(detector  signal)  per  cycle.   However,  16  y/sec  were 
needed  to  analyze  an L  record  each  incoming  signal.   A 
correction  for  this  deadtime  was  included  in  the  computer 
program  used  to  process  the  data  (see  Section  4. 3).   More 
de bailed  information  on  the  211  unit  and  the  multi-channel 
analyser  is  give?n  in  Reference  (4-). 

3.5.4-   The  Monitoring  and  Normalizing  System 

To  insure  that  the  pulsing  system  was  operating  under 
the  best  possible  conditions  at  all  times  a  beam  pulse 
monitor  v/as  provided. 

The  beam  pulse  monitor  was  necessary  due  to  the  pro- 
blem of  focusing  the  proton  beam  from  the  accelerator. 
Because  of  various  changes  and  drifts  in  the  accelerator, 
the  setting  for  optimum  focusing  changed  from  day  to  day, 
and  had  to  be  constantly  adjusted.   In  pulsed  operation, 
the  only  practical  method  of  focusing  was  to  use  an  electri- 
cal signal  from  the  target  itself  and  to  observe  the  signal 
with  a  fast  oscilloscope.   The  system  employed  is  sketched 
in  Figure  3.7.   Notice  that  the  target  was  insulated  from 
the  beam  tube  with  a  glass  ring  which  enables  the  signal 
to  be  picked  off  the  target.   The  monitoring  system  was 
especially  susceptible  to  R.F.  pick-up  from  the  pulser. 
All  cables  of  the  system  were  shielded  against  pick-up^ 
but  the 'problem  was  never  completely  eliminated.   A  typical 
target  puis")  can  be  seen  in  Figure  4.11  of  Vaughan's  thesis. 
The  electronic  pulse  being  delivered  to  the  deflection  platen 
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could  also  1    '  litored.'  This  system,  is  also  shown  in 
Figure  3»7« 

A  neutron  monitor  was  used  as  a  means  of  normalization 
for  comparing  different  runs  and  subtracting  background. 
A  detailed  discussion  of  normalization  is  given  in  Chapter 
5  of  Reference  (1).   As  in  previous  experiments,  this  author 
concluded  that  total  neutrons  produced  during  an  experiment 
would  be  the  most  reliable  quantity  to  use  of  normalization. 
A  boron  triflouride  counter  was  placed  between  the  graphite. 
stringers  of  the  assembly  about  13  inches  from  the  target 
and  on  the  side  opposite  the  absorber  foil  and  detector. 
Thus,  with  the  output  'of  the  BF^.  tube  going  to  a  scaler, 
the  time  integrated  neutron • flux  as  seen  by  the  detector 
could  be  reco, led. 
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CHAPTER  4 

Experimental  Results  and  Discussion 

•4.1 The  Use  of  Resonance  Absorbers  as  Neutron  Density 

Indicators 

Resonance  absorbers  were  used  in  this  experiment  to 
sample  neutron  population' at  the  resonance  energy.   The 
most  suitable  materials  for  this  purpose  are  those  .which 
exhibit  a  strong  isolated  absorption  resonance.   The 
higher  and  narrower  the  resonance,  the  better  it  approxi- 
mates a  delta  'function.   Table  II  of  Reference  (1)  gives 
a  list  of  absorbers  which  exhibit  strong  isolated  reson- 
ances.  The  energy  and  peak  cross  section  of  the  reson- 
ance" is  also  given. 

Another  important  consideration  in  choosing  a  suit-  : 
able  absorber  was  its  prompt  neutron  capture  gamma-ray 
•spectrum.   This  did  not  prove  to  be  a  problem  since  the 
de-excitation  gamma  rays  of  most  materials  are  spread 
over  a  wide  range  of  energy.   Reference  (1)  gives  capture 
gamma-ray  spectra  for  the  three  absorbers  used  in  these 
experiments:   cadmium,  indium,  and  gold.   The  capture  gamma 
spectra  of  these  materials  extends  from  the  KeV  region 
to  energies  exceeding  7  MeV. 

Indium  and  gold  were  chosen  for  their  strong  isolated. 
resonances  at  1.46  eV.  and  5*0  eV  respectively.   Cadmium,- 
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broad  resonance  that  reaches  a  maximum  at 
: about  C.2  eV,  was  also  used,  although  it  was  expected  that 
the  time-  resolution  for  these  measurements  v/ould  not  be 
as  good  as  that  for  indium  and  gold. 

4.2  Performance  of  the  Experiment 

In  order  to  insure  a  constant  geometric  relationship 
between  the  foil  and  detector,  a  foil  holder  was  used. 
The  fori  holder  consisted  of  a  one  inch  thick  slab  of  gra- 
phite, which  was^sized  to  fit  into"  the  channel  of  the  assem- 
bly and  was  bonded  to  the  end  of  the  detector  with  epoxy 
cement.   The  foil  was  taped  to  this  slab,  and  the  whole 
detector-foil  array  was  inserted  into  the  graphite  assem- 
bly through  the  channel  (see  Figure  4.1).        n 

To  measure  the  spatial  dependence  of  the  slowing  down 
time,  it  was  necessary  to  be  able  to  move  the  detector- 
foil  array  along  the  channel.   Graphite  slabs  of  varying 
thickness,  which  v/ere  inserted  in  the  channel  between  the 
source  and  foil,  assured  that  there  would  always  be  solid 
moderating  material  in  front  of  the  foil. 

Before  performing  the  actual  experiments,  short-  runs  ■  ■, 
were  made  while  varying  the  discriminator  of  the  Ortec 
413  Strobbed  Discriminator.   For  each  foil,  the  discrimin- 
ator setting  which  resulted  in  the  best  foil  counts  to 
background  .counts  ratio  was  determined. 

A  typical  set  of  experimental  data  is  shown  in  Figure 
4.2(a).  .These  results  are  for  a  7  cm  x  7  cm  x  0.05  cm  gold 
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foil  positioned  28  cm  away  from  the  target.   The  beam 
pulse  width  was  200  usee,, and  the  time  interval  between 
bursts  was  ^>Q0   usee.   The  channel  width  was  2^0  nsec. 
A  time  distribution  with  the  gold  absorber  in  place  was 
taken,  and  a  distribution  without  the  absorber  was  taken 
immediately  afterwards.   The  BF-,  counts  for  the  gold  foil 
was  251,111  and  for  the  background  was  259,992.   Notice 
the  very  prominent  peak  around  channel  20.   This  peak  is 
called  t.he  prompt  peak  and  occurs  when  the  burst  of  pro- 
tons hits  the  lithium  target.   This  peak  was  used  as  an 
accurate  zero  time  reference  in  the  experiment.   The  right 
half  of  the  graph  (notice  the  change  of  scale)  shows  the 
region  in  which  the  most  probable  slowing  down  time  occurred 
No  peak  is  discernible  from  the  raw  data. 

4.3   Normalization,  Background  Subtraction,  Dead  Time 

Correction  and  Smoothing 

As  mentioned  in  Section  3«5«4,  the  time  integrated 
neutron  flux,  as  derived  from  a  BF-,  tube,  was  used  for 
normalization.   The  number  of  31?   counts  registered  by 
the  scaler  was  recorded  for  each  experiment.   A  factor 
obtained  by  taking  the  ratio  of  background-run  BF^  counts 
to  foil-run  BF-,  counts  was  aoolied  to  the  foil  data. 

In  order  to  obtain  the  net  effect  of  neutron  cap- 
tures in  the  foil,  the  counts  originating  from  all  other 
events  must  be  subtracted  when  the  time  distribution  is 
taken.   These  sources  of  parasitic  gamma  rays  include: 
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1)  Inelastic  scattering  by  the  target,  the  absorber 
foil,  and  the  detector. 

2)  ■  Absorption  in  materials  other  than  the  absorber 

* 

foil. 

3)  Reactions  coincident  with  the  neutron  burst. 

i 

4)  Overlapping  events  from  preceding  bursts. 

5)  Steady  state  background. 

All  these  items  were  taken  into  account  by  taking 
a  time  distribution  without  the  absorber  foil  in  place, 
and  then  subtracting  it  from  the  distribution  taken  with 
the  foil. 

As  mentioned  in  Section  3»5»3i  there  was  a  16  "usee 
dead  time  associated  with  each  count  registered  in  the 
multi-channel  analyzer.   In  Appendix  C  of  Reference  (2) , 
a  theoretical  discussion  of  the  dead  time  correction  for 
the  TMC  Model  211  Unit  is  presented.   The  correction  fac- 
tors derived  in  that  reference  were  used  in  this  thesis. 

Computer  programs  which  normalized,  subtracted  the 
background,  and  applied  dead  time  corrections  to  the  raw 
data  are  listed  in  Appendix  A. 

Figure  4.2(b)  shows  the  time  distribution  resulting 
from  processing  the  raw  data  in  Figure  4.2(a)'.   It  is 
obvious  from  this  graph  that  there  is  a  maximum  in  the 
time  distribution  occurring  roughly  at  channel  number  60. 
It  is  also  obvious,  however,  that  the  fluctuations  in  the 
curve  prevent  an  accurate  determination  of  the  maximum 
point.   To  reduce  this  uncertainty,  a  computer  program 
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at  M.I.T.  (3)'  which  applied  a  smoothing  tech- 
nique to  the  corrected  data  was  used.   Effectively,  the 
program -fits  a  Fouier  series  curve  to  the  data  points. 
Smoothing  is  accomplished  by  choosing  not  to  fit  the  ser- 
les  toj&igh  frequency  waves,  which  correspond  here  to  sta-  . 
tisticfu  fluctuations.   The  frequency . cutoff ,  above  which 
all  hiwher  frequencies  will  be  thrown  away,  is  under  the 
contriJr  of  the  programmer.   Thus,  any  desired  degree  of 
smooti"»rg  can  be  effected  to  a  set  of  data  points.   In 
this  experiment,  all  frequencies  less  than  four  channels 
wide  were  dropped.   Figure  4.2(c)  is  the  smoothed  version 
of  Figure  4.2(b).   Using  this  smoothed  distribution^  it 
was  estimated  that  the  maximum  point  of  the  time  dis- 
tribution occurred  in  channel  62.5+1.   Since  the  prompt 
peak  is  in  channel  20,  the  slowing  down  time  is:  (43.5+1) 
x  (.25  ]i sec/channel)  =  10. 87+. 25  usee.   The  smoothing 
program  is  listed  in  Appendix  A. 

4.4  Space  and  Energy  Dependent  Slowing  Down  Time  Mea- 
surements 

If  various  absorbers  with  resonances  located  at 
different  energies  are  used  and  if  the  distance  of  these 
absorbers  from  the  source  is  varied,  the  time  dependent 
neutron  distribution  can  be  measured  as  a  function  of 
space  arid  neutron  energy. 

Figures  4.3  to  4.6  represent  time  distributions  to 
the  cadmium  resonance',  Figures  4.7  to  4.11  represent  time 
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distribution to  the  indium  resonance^  and  Figures  4.13 
to  4.16  represent  time  distributions  to  the  gold  reson- 
ance.  In  each  Figure,  (a)  represents  the  distribution 
before  smoothing  and  (b)  after  smoothing.   Notice  that 
the  scale  for  the  cadmium  curves  is  500  nsec/channel , 
while  it  is  2f?0  nsec/channel  for  the  indium  and  gold 
curves. 

The-  results. of  this  experiment,  together  with  their 
probable  errors,  are  graphically  presented  in. Figures 
4.17  to  4.19.   The  experimental  uncertainty  associated 
with  each  point  reflects  the  broadness  in  the  time  dis- 
tribution- around  the  maximum  value.   As  expected,  the 
peaks  of  the  cadmium  distribution  were  much  broader  than 
those  of  gold  or  indium'.   The  solid  line  in  the  gold  and 
indium  figures  has  been  constructed  from  theoretical  time 
distributions  (see  Section  2.3),  based  on  diffusion  theory 
(4). 

The  experimental  results  for.  the  slowing  down  time 
to  5«0  eV  appear  to  be  in  good  agreement  with  the  theo- 
retical calculations.   The  agreement  with  theory  for 
1.46  eV  is  not  as  good.   The  experimental  results  indi- 
cate that  there  is  little,  if  any,  spatial  dependence  of 
the  slowing  down  time  for  distances  closer  than  about 
30. cm  from  the  source.   For  distances  greater  than  30  cm, 

the  experimental  results  show  a  sharp  rise  in  t  with  dis- 

m 

fcance,  indicating  a  strong  spatial  dependence.   Theore- 
tical calculations  do  not  predict  as  sharp  a  rise  of 
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t  as  shown  by  the  experiments. 

This  disagreement  "between  theory  and  experiment  was 
probably  due  to  the  breakdown  of  the  infinite  medium 
boundary  conditions,  since  the  foil  was  only  9  cm  from 
the  outer  edge  of  the  assembly  at  the  last  spatial  point. 
The  absence  of  moderating  material  surrounding  the  foil 
meant  that  the  neutrons  took  a  longir  time  to  slow  down 
to  a  given  energy.   More  measurements  need  to  be  taken 
at  greater  distances  to  resolve  this  question. 

A  few  measurements  were  made  with  the  detector  separ- 
ated from  the  foil  by  about  10  cm  of  graphite  in  order 
to  determine  if  the  perturbation  caused  by  the  detector 
had  any  effect  on  the  slowing  down  time.   These ^.results 
are  marked  by  triangles  in  the  Figures.   It  can  be  seen 
that,  the  gold  and  cadmium  results  agree,  within  the  limits 
of  experimental  uncertaintity,  with  previous  measurements. 
The  indium  point,  however,  falls  far- below  the  original 
measurements.   The  author  can  give  no  explanation  for  the 
contradiction  of  this"  point  with  either  the  original  results 
for  indium  or  the  close  agreement  between  the  first  and 
second  measurements  for  the  gold  and  cadmium.   Further 
measurements,  with  the  detector  and  foil  separated,  need 
to  be  performed  to  clear  up  this  contradiction. 

4-.  5   Conclusions 

■  There  was  some  agreement  between  the  spatially  depen- 
dent slowing  down  time  measurements  made  in  these  experi- 
ments and  theoretical  calculations  using  diffusion  theory. 
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No  measurable  spatial  dependence  of  the  slowing  down  time 
was  observed  for  distances  less  than  30  cm.   Spatial  depen- 
dence was  observed  for  distances  greater  than  30  cm,  although 
this  dependence  was  much  more  promounced  in  the  experiments 
than  predicted  by  slowing  down  theory.   More  space  depen- 
dent measurements  at  distances  greater  than  4-4-  cm  should 
be  made  in  a  larger  graphite  assembly  to  determine  the 
true  spatial  dependence.   Measurements  to  determine  'if 
the  perturbation  of  the  detector  had  any  effect  on  the 
slowing  down  time  were  inconclusive  and  should  be  repeated. 
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Appendix  A.   Computer  Program  Listings 

A.l 
C      FORMALIZATION,  DEAD  TIME  CORRECTION,  AND  BACKGROUND 
C      SUBTRACTION  CODE  LISTING 
C 

DIMENSION  T8M(2  55),CHAN(255) , BACK { 255 ), FOIL! 255) f ANORM( 255 ) , 
1C8ACK(255),CF0IL(255),X(255) , Y ( 2  55 ) , TBAC 20  6 ) , ACHAN t 206) ,B{800) 
COMMON  TBM, CHAN, WIDTH 

EQUIVALENCE(TBM(50JtTBAtl)!t ( CHAN ( 50 ), ACHAN ( 1 ) ) 
READ  2  00, WIDTH 
200  FORMAT(F6.2) 

READ  99,8C0UNT,FCQUNT 
99  FORMAT  <2F 10.0)  / 

NUM=0  / 

103  NUM=NUM+1 

READ  100,  (BACKt I ) ,1  =  1,255; 

100  FORMAT  (3X,10F7.0) 

READ  101,  (FUIKI  ), 1=1, 255) 

101  FORMAT  (3X,10F7.0) 

READ  98,3BUR,ABUR  i 

98  FORMAT  12F10.0) 

PRINT  701,BC0UNT,FC0UNT  i 

701  FORMAT  (3X,2F10.0> 

PRINT  100,  (BACKd  )  ,1  =  1,255) 
PRINT  101»IFOILU  ), 1  =  1, 255) 
PRINT    702,BCUR,ABUR 

702  FORMAT     13X,2F1Q.Q) 
L  =  -l 

420    L  =  LH 

IF    ID419,  419,417 
419    DO    418    1=1*255'  \  \ 

418    X(  I  )=BACK(I) 

BUR=BBUR 

GO  TO  415 
417  DO  416  1=1,255 
416  X(I )=FOIL(I) 

BUR  =  ABUR  :.\,  ! 

415  0UM  =  1.6. /WIDTH  •  ' 

N  =  DUM 

YM)=BUR*LQGF(BUR/(BUR-X(1)  )  )  \ 

SUM=0.  , 

DO  610   I=2,N 

SUM=SUM+X( 1-1)  •         i 


-30- 


610 


615 

414 
413 

412 

All 


906 

306 

50 


Y{  I  ) 
SUM  = 

NM 
DQ  6 
Y(I) 
K  =  I- 
SUM  = 
IF  { 
DO  4 
CBAC 

GO 
DQ  4 
CFOI 
i  =  BC 
DO  9 
AN  OR 
TBM( 
CHAN 
CONT 

prin 

FORM 


=  BUR 
15  I 

=our< 
N 

SUM  + 

LJ41 
13  I 
K(  I  ) 
TO  4 
11  I 
L(  I  ) 
OUNT 
06  I 

Mm 

I)=A 
(I)  = 
INUE 
T  30 

AT  I 


»L0GF(  { BUR-SUM) /(BUR-SUM-XI  I))  ) 
X  t  N ) 

=MM,255 

»L0GF( ( BUR-SUM) /(BUR-SUM-X( I )> ) 

X( I )-X(KJ 
4,414,412 

=1,255 

=  YU) 

20 

=1,255 

=  Y(I) 

/FCOUNT 

=1,255 

=  Z»CFQ[L(I  ) 

NORM (I  J-CtJACK(I) 

floa'tfu) 

6, ITBM(I), 1=1,255) 

1H  ,10F10.1) 


CALCOMP  ROUTINES 

CALL  PL0TS1(B,800) 

CALL  SYM3L5(3.0,9.0,D.2,20H  ID  NO.  =  M3994-5788, 0. , 20 ) 

CALL  NUMBRl{4.0,8.5,0.2,NUMf0«.t-U 

CALL  PICrUR[8.5,4.2,15H  CHANNEL  NUMBER, 

15,19H  COUNTS  PER  CHANNEL , 19, ACHAN, TBA, 206, 0. , 0 ) 

END  FILE  10 


700 


READ  99, 3C0UNT, FCOUNT 
IF  (BCOUNT)700,700,103 
CALL  EXIT 

END 
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A. 2 

C 
C 

c 


DATA 
TIME 


SMO 
COR 


OTHING  CODE 
RECTI  ON,  AND 


LISTING  INCLUDING  NORMALIZATION,  DEAD 
BACKGROUND  SUBTRACTION 


200 
102 

201 

205 
203 
99 
103 

100 

101 

98 

701 

702 
420 

419 

418 


417 
416 


OlMENSIO 
1XC0S4137 
28ETA(75) 

3CBACM25 
COMMON  T 
UIYALE 
READ  200 
FORMAT  { 
READ  102 
FORMAT  ( 
PRINT  20 
FORMAT  ( 

15H  CRT=, 
PRINT  20 
FORMAT  ( 
PRINT  20 
FORMAT  ( 
READ  99, 
FORMAT  { 
NUM-0 
NUM=NUM+ 
READ  100 
FORMAT  { 
P.E\D  101 
FORMAT  ( 
READ  98, 
FORMAT  ( 
PRINT  70 
FORMAT  ( 
PRINT  10 
PRINT  10 
PRINT  70 
FORMAT  { 
L=-l 
L=L+1 
IF  U)41 
DO  418  I 
X(I )=BAC 
t=BBUR 
GO  TO  41 
DO  416  I 
I )=FOI 
BUR^ABUR 


N  TB 

J,TB 

,GAM 

5)  ,C 

8N,X 

NCE( 

tlW, 

I  4 ,  E 

,  IWT 

12F6 

1 ,  I  W 

1H1, 

F6.0 

5 

1H  , 

3 ,  I  W 

1M  , 

BCQU 

2F10 


N(37 
MI  25 

MA  (7 

FOIL 

S IN , 

tb:<{ 

OMEG 

10.4 

(I), 

.3) 

,OME 

9H  N 

tiX, 


GAH, 
-OME 
8H  F 


CRT, FACTR, WIDTH 

GA=,I4,1X,9H  D-QMEGA=,E1Q.4,1X, 

ACT0R=,F6.3,7H  WIDTH=,F6.2) 


20X,16H 
Tf 11,1*1 

5E15.7) 

NT,FCOUN 
.0) 


WEIGHT 
>  IW> 


FUNCTION) 


, [BACKCI > t 1= 

3X,10F7.0) 

,{FOILU) ,1= 

3X»10F7.0J 

BBURf ABUR 

2F10-0) 

1,BC0UNT,FC0 

3X,2F10.0) 

0, (BACMI) ,1 

1,  (FOIL*  I  )  ,1 

2,B8UR,ABUR 

3X,2F10.0) 


9,419,417 
=1,255 
K[  I) 

5 

=1,255 

LCI) 


1,255} 
1,255) 

UMT 

=1,255) 
=1,255) 


■  .* 
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415 


610 


615 

4L4 
413 

412 
411 


905 


902 
901 


903 


906 


DUM  = 
N  =  DU 
Yd) 

DO  6 
SUM= 

y  •:  i  ) 

SUM  = 
NM 

00  6 

Yd ) 

K=I- 

SUM  = 
IF  { 

on  4 

CBAC 

GO 
DO  4 
CFOI 
DO  9 
YKAK 
XCOS 
CONT 
DO  9 
DO  9 
XKAK 
SEIG 
CONT 
CONT 
DO  9 

ALPH 

1C0SF 
BETA 

1SEIG 
GAMM 
CONT 
DQ  9 
Z  =  3C 

TBM( 
CONT 
DO  1 


16 

'I 

=  B 
0, 
10 

su 

—  p 

su 

=N 

15 

=  B 

N 

SU 

L) 

13 

K( 

TO 

11 

L( 

05 

U4 

41 

IN 

01 

02 

U= 

EN 

IN 

Ih 

03 

A{ 

A{ 

(T 

(J 

EN 

' 

IN 

06 

OU 

Hi 

I) 

IN 


./WIDTH 

UR*.LOGFIBUR/(BUR~Xm ) ) 

1  =  2,  N 
M  +  X( I- I) 

LOGF{ (BUR-SUM) /{BUR-SUM-Xt I ) ) ) 
H*-X(NJ 

«■] 

I=NM,255 
UR*LOGFl (BUR-SUM)/(BUR-SUM-XU ) ) ) 

M+X{ r )-XlK) 
414,414,412 

1=1,255 
I )=Y( I ) 

420 

1=1,255 
I)=Y(I) 

1=1,37 
=  C"EGAH*FL0ATF(IW-1)»FL0ATFU) 
I )=C0SF(YKAKU4) 
UE 

1=1,37 

J=1,IW 
OMEGAH*FLOATF { I ) *FLOATF{ J-l ) 
(I, J)=SINF(XKAKU) 
UE 
UE 

J=1,IW 
J ) =ONEGAH*FLOATFl J-l ) 

JI=THETA{ J)**l-3)*ITHETA(J)*«2*THETA{ J ) »SEIGEN C 1 , J 
HETA(J))-2.*SEIGEN(1.J)«*2) 

)=2.*THETAIJ)»*(-3)*(THETA( J)*(1.*C0SF(THETA( J) )** 
U,J)*COSF{THETAU))  ) 

J)=4.*THETA(J) »*t-3)*(SElGEN<lf J)-THETA( J)*C0SFITH 
UE 

1=1,255 
NT/FCOUNT  '     . 

n-z»CFoiLm 

=ANORM(I)-CBACK{I J 

UE 

1=1,27 


i 
2)-2.* 

ETAU) ) ) 
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90  \ 


511 


517 


DO  904 
YKAKU1 
YKAKU 
XCQSK 
XC0S21 
CO NT  IN 
5 1 1 
YA=ALP 

YOG  AM 
SXG4=0 

sic 

DO    1    N 

SIG3=S 

DO    2    N 

SIG4=S 

XA=YA* 

XSINU 

CONTIN 

DO    517 

YA=ALP 

YB=BET 

YC=GAH 

SIG^O 

SIG3=0 

NY=IW- 

DO    3    N 

SIG3=S 

NY=IW- 

DO    4    N 

SIG4=S 

XAOME 

IN=IX 

TSK.i  IN 

CHANt I 

CONTIN 

DO    914 

YKAKU1 

YKAKU2 

XCOSK 

XCOS2( 


J=1,IW 
=  0,MEGAH*FLDATF(J-1) 
=27.»YKAKU1 
J)=COSF(YKAKUI) 
J)=C0SF{YKAKU2) 
UE 

1=1, IW 

HAin 

A(  I  ) 
MAI  I  ) 

«5MTBNm*SEIGENU,I)+TBNI27)*SEIGEN{27,I>) 

N=3,25,2 

IG3+TBN(NN}»SE!GENINN,I) 

N=2,26,2 

IG4  +  TBN(NN)»SElGEN(N.M,n    ' 

(T8N{1)»XC0S1(  I)-TBN(271*XC0S2(  I  )  )  *YB»S  I  G3*YC*S  I  G4 

J=XA*WT(I) 

UE  li 

IX=1,16 
HA(IX*1) 
A( IX+1) 
MAUX+1) 

.5*IXSIN(1) *SEIGEN(IX,  U  +  XSINl I W ) *SE I GENt IX, I W  )  ) 

2 

N=3,NY,2 

IG3+XSIN(NN)*SEIGENCIX,NN) 

1 

N=2,NY,2 

IG4+XSIN(NN)*SEIGENl IX,NN) 

GAH»iYAMXSlNt  l}-XSINUVn*XC0S4(IX)  )  +YB*SIG3+YC*SI  G4  ) 

)=Q.63662*XA 
N)=FLQATFC IN) 
UE 

J  =  1 ,  I  W 
=OMEGAH*FLOATF(J-l) 
=37„«YKAKU1 
J)=C0SF(YK4KU1) 
J)=C0SF(YKAKU2) 
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914 


525 


521 


527 
530 


535 


c  o  n  r  i  ^ 

DO  530 
DO  525 
K-16«( 
TBN1  I  ) 

DO 

YA=ALP 
YB  = 

YC=GAM 
SiG4=0 

3  =  0 
DO  5  N 
SIG3=S 
DO  6  N 
S1G^=S 
XA=YA* 
XSIN{ I 
CONTtN 
DC  5  27 
YA=ALP 
YB=BET 
YC=GAM 
SIG4=0 
SIG3=0 
NY=IW- 
DQ  7  N 
SIG3=S 
NY=IW- 
DO  3  N 
SIG4=S 
XA-QME 
IN=16» 
TBxt  IM 
CHAN t I 
CONTIM 

TIN 
DO  535 
K=230+ 
I  ::![[) 
DQ  924 
YKAKU1 
YKAKU2 
XCO.SK 
XC0S2( 


UE 

J  =  2 

1  =  1 

J-l) 

=TBM 

1  =  1 

HA  (I 

All) 

MAI  I 

.5  I 
N=3, 

IG3  + 
N=2, 
IS'n- 
tTBN 

)=XA 
UE 

IX  = 
HA(  I 
At  IX 
MA  {I 


,15 
.37 

-104-1 

(K) 

,  IW 
) 

) 

TBN« 1)*SEIGEN(1, I >+TBN( 37) *SEI GEN137, I) > 

35,2 

TBN(NN)*SEIGEN(NN,I> 

36,2 

TB.vUN.\)*SEIGEN(NNtI) 

{1)*XC0S1(I  )-TBN{37)*XC0S2CI)  ) *YB«S I G3+YC*SI G4 

*WT(I) 

11,26 

X  +  1J  U 

+  1) 

x+n 


.5MXSINU)*SEIGEN{IX,l)«-XSIN[IW)»SEIGEN(IXf IW) ) 

2 

N  = 

IG 
1 

M  = 


3, NY, 2 

;3+XSIN{NN)»SEIGEN(IX,NN) 


GA 
(J 
)  = 

H) 
UE 
UE 

I 
I 
=  T 

J 
=  0 
=  2 
J) 
J) 


2, NY, 2 

4  +  XSlN(NN)»SEIGENnx,NN) 

H* (YA*(XS1N(1)-XSIN(IW)»XC0S4MX) ) +YB»SIG3+YC*SI G4 ) 

-D-HX-10 

0.63662*XA 

=FLQATF(IN) 


=  1,25 

BM{K) 

- 1 ,  I  W 

KEGAH*FL0ATF(J-1) 
5.»YKAKU1 
-COSF(YKAKUl) 

=COSF{ YKAKU2) 
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92'+      CONTINUE 

DO    531     1  =  1,  I W 
YA=ALPHA(I) 
YB-BETAU) 
YC  =  GAMMAU) 
SIG'^0. 

SiG3  =  0.5*{TBN(l)  »SEIGENU,I)«-TBN(25)*SEIGEN{25,I  5  ) 
DO    9    NN=3,23,2 
9      SIG3=SIG3+TBN{NN)«SEIGEN(NN,n 
DG    10    NN=2,2Vt2 
10      SIG4=SIG4+TBN(NN)*SEIGEN(NN, I) 

XA  =  YAMTBNC 1)«XC0S1U ) -TBN { 2 5 ) *XC0S2 ( I) )  «-YB*S  IG3*YC*SI  G4 
XSINU>=XA*WT<I) 
531      CONTINUE 

DO    537     IX=ll,25 
YA&ALPHA1 IX+1) 

A  IX+1) 
YC=GAMMA(IX+1)  j 

SIG4=0. 

SiG3  =  0.5*{XSINll)*SEIGENMXf  1)*XSIN(  IW)*SEIGEN(  IX,  IW)  ) 
NY*lW-2 

DO    11    NN=3,NY,2 
'      11       SIG3  =  5IG3  +  XSIN{NN)*SEIGENUXfNN) 
NYsIW-1 

DO    12    NN=2tNY,2  ! 

12      SIG^=SIG^+XSIN(NN)*SEIGEN(IX,NN)  j 

XA  =  0MEGAH*IYAMXSIN{1>-XSIN(IW)»XC0S*IIX) ) +  YB*S I G3*YC*S I G4 ) 
IN=23Q+IX  I 

TBK( In)=0.63662*XA 
CHANUN)*FLOATFtIN) 
53  7       CONTINUE 

PRINT  307 
307   FORMAT  (1H1,1SH  SMOOTHED  FUNCTION) 

PRINT  306, ITBXil ), 1=1, 255) 
306   FORMAT  UH  ,  10F10.1) 
C 

C      CALCOM?  ROUTINES 
C 

50   CALL  PLOTS! (Br  300) 

CALL  NUMBRH<r.0,8.5,0-2tNUM,0.  ,-1) 

LL  SYMBL5I3.0,9.0,0.2»20H  ID  NO.  =  M3994-5788, 0. , 20 ) 
CALL  PICTUR(8.5,4.2, 15H  CHANNEL  NUMBER, 
115.28H  SMOOTHED  COUNTS  PER  CHANNEL , 28 , ACHAN, TBA, 21 6, 0. , 0 ) 
END  FILE  10 
C 

READ  99tBC0UNT,FC0UNT 
IF  tBCOUNT) 700,700,103 
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700  CALL  EXIT 

END 
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